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Abstract 

The electron transfer between an F~ ion and Cm(111) and Ag{lll) surfaces is studied by 
the wave packet propagation method in order to determine specifics of the charge transfer 
interaction between the negative ion and the metal surface due to the projected band gap. A 
new modeling of the F~ ion is developed that allows one to take into account the six quasi- 
equivalent electrons of F^ which are a priori active in the charge transfer process. The 
new model invokes methods of constrained quantum dynamics. The six-electron problem is 
transformed to two one-electron problems linked via a constraint. The projection method 
is used to develop a wave packet propagation subject to the modeling constraint. The 
characteristics (energy and width) of the ion F~ ion level interacting with the two surfaces 
are determined and discussed in connection with the surface projected band gap. 
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I. INTRODUCTION 



When an atom or a molecule is close to a surface of a solid, its electrons interact with those 
of the solid, leading to the possibility of an electron transfer between the atom (molecule) 
and the solid. This charge transfer process plays a very important role in a variety of 
different situations. In particular, it often occurs as an intermediate step in reactions at sur- 
faces (desorption, fragmentation of adsorbates, chemical reactions, quenching of metastable 
species, etc.) The one-electron transfer between energetically degenerate electronic 

levels of the atom and the solid is called the Resonant Charge Transfer (RCT) process. 
It is usually considered as the most efficient one among various possible charge transfer 
interactions. Since a few years the development of accurate theoretical approaches to the 
RCT in the case of free electron metal surfaces [p|-[TT|| has led to a successful description 
of a one-electron transfer in the interaction of ions (atoms) with such surfaces |1T^-[14|. All 
these approaches are based on the description of single electron being transferred between 
the atom and the surface. As an example, one can mention the neutralisation of an alkali 
positive ion by RCT, even if the possibility of capturing the electron in two different spin 
states somewhat alters the one-electron picture |l^,|T6|. However, many atoms or molecules 
contain more than one electron which can possibly participate in the charge transfer pro- 
cess, especially, when the active electrons occupy the same energy level. The latter has to 
be taken into account in any quantitative approach to the RCT 

For instance, a free fluorine negative ion can be described as a closed 2jP outer electronic 
shell with six equivalent electrons, and all of them can participate in the RCT when the ion 
is close to a metal surface. The effects of each electron cannot be simply added up to get 
a total effect. Indeed, consider the loss of an electron by the ion. Any of the six electrons 
can be the one that is detached. However, the detachment of one electron a 'priori precludes 
the detachment of another one. Clearly, the loss of a second electron would be a completely 
different process because it corresponds to a formation of a positive fluorine ion and, hence, 
has a different energetics. Such correlations between the outer-shell electrons of the F~ 
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ion must be accounted for in any description of the RCT. In this work we show how such 
correlations can be modeled by two one-electron systems linked via a constraint. 

In classical dynamics, constraints appear as some algebraic relations between generalized 
coordinates of a system and their time derivatives (velocities), which are to hold for any mo- 
ment of time. They are widely used to model, e.g., effects of an environment on a system 
in question, or to develop theories with local symmetries (gauge theories) such as electro- 
dynamics, general relativity and Yang-Mills theories. Classical constrained dynamics has 
been well studied in classical works by d'Alambert, Lagrange, Holder and Gauss (see, e.g., 
Ref. [|l^ and references therein). In quantum mechanics constraints have been analyzed by 



Dirac |]T9[, mainly because of the need of quantizing fundamental physical theories with local 



symmetries. The question addressed by Dirac was the following: Given a classical system 
with constraints, construct a quantum theory which satisfies the correspondence principle. A 
time evolution of a quantum system is described by the evolution operator. Its kernel in 
the coordinate representation is called a quantum mechanical propagator. For constrained 
systems it is usually obtained by a reduction of the Feynman phase space path integral onto 
the physical phase space, as proposed by Faddeev [Q. There are many subtleties associated 
with quantization of constrained systems (see for a review, e.g., [pi|]). 

Here we propose a quantum constrained system which can be used to model the charge 
transfer interaction between an F~ ion and a metal surface. The basic idea is to transform 
the original six-electron problem to two one-electron problems each of which describes a 
possible decay channel of the ion. As a consequence of the quantum equivalence of the 
outer-shell electrons, the effective one-electron systems turn out to be linked only by a 
constraint. In other words, the entire interaction between the two systems occurs through a 
constraint rather than via a local potential. One can also say that such a kinematic coupling 
of the decay channels is induced by the symmerty of the original six-electron problem. 

The constraint has a remarkable feature: It does not have a classical limit, meaning 
that its effects disappear in the formal classical limit h 0. That is, the dynamical effect 
modeled by such a constraint is essentially quantum and cannot be modeled by any classical 



potential force. In this regard, the canonical quantization scheme for constrained systems 
due to Dirac is not applicable here. We develop below a novel computational scheme for 
a quantum evolution subject to such purely quantum constraints. The scheme is based 



on the projection formalism introduced earlier by one of us I^Sj within the framework of 
gauge theories. Our approach also comprises a novel method to account for the intra-atomic 
correlations within a one-electron description of the charge transfer interaction between an 
F~ ion and the C'u(lll) and Ag{lll) surfaces. We demonstrate that the approach turns 
out to be very efficient in the wave packet propagation studies. 

The choice of the F~ / Cu{lll) and / Ag{lll) systems is motivated by several reasons. 
The interaction of halogen negative ions with a free-electron metal surface has already been 
studied theoretically within the Coupled Angular Mode approach which lead to a quite 
satisfactory description of the halogen negative ion formation in scattering halogen atoms 
by various metal surfaces [|12|,|13|. Recently, on the example of the Cm(111) surface it has 
been shown that the projected band gap of the (111) surfaces of noble metals strongly affects 



the RCT p5|P^. The point is that in a certain energy range electrons cannot propagate 



along the normal to the surface (L-band gap in the < 111 > direction On the other 

hand, the RCT process corresponds to an electron tunneling between the atom and the 
surface which is favored along the surface normal. In addition, it was shown that the 2D 
electronic continuum of the surface state plays a significant role, often dominating the RCT 
process. The band gap has several consequences. First, there exist very long lived states in 
the alkali-CM(lll) systems p4| , p6| , p7| . Second, there is a parallel velocity dependence of the 



probability of an electron capture by the atom from Cm(111) surfaces in grazing scattering 



experiments that is characteristic of a 2D electronic continuum [^. Third, there exits an 
avoided crossing between the energy level of the projectile and the bottom of the surface 
state continuum |23]. These effects have been found to depend strongly on the interaction 



time p3| , p9| , p0| : They only appear for long interaction times, i.e., for slow collisions. 

The free F~ ion energy level is slightly above the bottom of the surface state continuum 
of the Cm(111) and Ag{lll) surfaces, and therefore, because of the image charge attraction. 



it could cross the bottom of the 2D continuum for a finite ion-surface distance. In the present 
work, we investigate the effects of the projected band gap and of the surface state continuum 
on the F~-metal charge transfer. In particular, we study the behavior of the system when 
the ion energy level is very close to the bottom of the 2D surface state continuum. A Wave 
Packet Propagation (WPP) approach to the RCT PJ23[] provides a quantitative description 



of dynamic and static aspects of the F -surface charge transfer. The latter allows one to 
analyze the dependence of the band structure effect upon an interaction time in the RCT. 

II. METHOD 

A. A negative ion F~ 

The free negative ion F~ is usually described in the Hartree-Fock approximation as a 
closed shell ion with the electronic configuration ls^2s^2p^. Its binding energy is 3.4 eV. In 
this approach the outer shell electrons are regarded as equivalent. Their wave function is 
given by the corresponding Slater determinant 

$ = ^ 2K2p^2p?2pf2p^i2/i , (1) 
V d! 

where 2p symbolizes a 2p orbital wave function, the superscripts a and (3 stand for the 
two possible spin directions, and the subscript indicates the magnetic quantum number m 
corresponding to the projection of the 2p angular momentum on the quantization axis. The 
Slater determinant (|l|) can be expanded into a sum of products of wave functions of an ionic 
core and an outer electron: 

^=\j:\FAA,ct^, = j:\G,\ct^,, (2) 

where the five-electron determinant | Fj\ describes a state of the fluorine atom F(^P). The 
product AjCpj corresponds to the 2p^'^ orbital, which has been factorised into a spin factor 
Aj and a spatial wave function 0j = (f)j{r). The wave function can also be used for 
an open shell description of the negative ion of the type 2p^2p'. In this case, the spatial 



wave function (pj singled out in corresponds to the outer 2p' orbital of the negative ion, 
whereas the core wave function \Gj\ is formed by the inner 2p orbitals. When analyzing the 
electron detachment process in the open shell description, the 2p' orbital is regarded as an 
active one. 

The representation (0) is particularly well suited for an analysis of the electron cap- 
ture/loss process between a fluorine ion F~ and a metal surface since it gives an expansion 
of the ion wave function over possible detachment channels. So, we retain the representation 
(I) to describe a loss or capture of an electron by the fluorine core 

^ = ElG',lV^. . (3) 
J 

Here ipj = V'i(^) is a wave function of an electron (captured or lost) in the j-channel. Ne- 
glecting the spin-orbit interactions and assuming the ion-surface interaction to be invariant 
under translations of the ion parallel to the surface, the system becomes invariant under 
the spin flip and rotations about the z-axis which is set to be perpendicular to the metal 
surface and passing through the ion center. Next, the z-axis is chosen as the quantization 
axis. Therefore the states with the 2;-component of the electron angular momentum ±1 are 
degenerate. As a consequence, only two electron wave functions are distinct in the repre- 
sentation d^). They correspond to the states with m = and |m| = 1. In what follows we 
assume that the charge transfer does not affect the neutral core wave function \Gj\. Only 
the outer electron is subject to the RCT dynamics. 

Any state of the system can be represented as a two dimensional isovector |\E') whose 
components are one-electron states corresponding to m = and |m| = 1 (upper and lower 
elements of the isovector, respectively). In cylindrical coordinates {z, p, (p), it can be written 
as 



(f |\&) 



(4) 



Note that, due to the symmetry of the problem, the (p dependence can be explicitly given. 
Hence, one can limit oneself to studying the z- and p-dependence of the electron wave packet. 



In our representation a free ion wave function is given by 




(5) 



where p{r) is the radial part of the free-ion orbital and Yi^iO, if) are the spherical harmonics. 
The different normalization factors of the isovector components are due to the fact that the 
free ion wave function contains twice as many states with |m| = 1 as with m = 0. The 
functions p{r)Yim{0 , ip) have a unit norm so that (\l'i?- |\E'ir-) = 1. 

Yet another remark is that the approach outlined above assumes that only an outer 
electron can be detached (the detachment occurs through the evolution of the wave functions 
ipo and ipi), i-e., it assumes an open shell description {2p^2p') of the ion F~ . A closed shell 
description {2p^) would correspond to an increase of the charge transfer interaction by the 
factor -\/6, and, hence, to an increase of the width by the factor 6 (see a discussion in [0). 
This appears to be better adapted for the halogen negative ion case. 

With all the above settings, a modeling of the charge transfer in the F~-metal system 
implies finding a one-electron Hamiltonian that governs a time evolution in the Hilbert space 
spanned by vectors (§). We take it in the following form 



The operators H, T, V^t, Vs and Hat are diagonal 2x2 matrices, in fact, we choose them 
to be proportional to a unit 2x2 matrix, with the diagonal elements denoted, respectively, 
if, T, Vat, Vs and i^at- Here, T is the electron kinetic energy operator, Kt the potential of 
the interaction between an electron and the neutral core, Vs the potential of the electron- 
surface interaction, i^at = T + K,t, and H = H^t + Vs. The wave functions p{r)Yim{0, ip) are 
eigenfunctions of the one-electron Hamiltonian iJat- 

The time evolution of the wave function (^) generated by the Hamiltonian (^ is nothing 
but the evolution of two independent one-electron wave packets. On the other hand, the two 
components of the isovector (^ cannot evolve independently. Indeed, in the free ion case, the 
two components cannot be arbitrarily chosen in order for the state (^ to describe a physical 
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H = T + Vat + Vs = Hat+Vs . 



(6) 



free ion. Actually, the radial components of the column elements appear to be proportional 
to each other with a specific factor (cf. (^). This relation between the two components 
comes from the expansion of the Slater determinant (|1|) which possesses a high symmetry 
being the symmetry of a quantum system of identical particles occupying the same energy 
level. The very same symmetry must be preserved in the expansion (|^) and, hence, upon 
a reduction of the six-electron description to our one-electron formalism. In other words, 
there must be a correlation between evolving components of the isovector corresponding to 
the electron states with m = and |m| = 1 thanks to the symmetry of the six-electron 
problem. Therefore physically admissible states in the Hilbert space spanned by isovectors 
(I) must be subject to some constraints required by the symmetry of the original six-electron 
problem. This will be a key point of our new approach to the RCT dynamics. 

To illustrate the necessity of constraints, consider the case when the detachment occurs 
in the m = channel. Then the bound part of the |m| = 1 channel must also disappear at 
the same time because there is only one ion F~ which contains both m = and |m| = 1 
components. In the Effective Range approach |^ used in the Coupled Angular Mode 
(CAM) method 0, this problem has been solved in the following way. One only considers 
the wave function of the system given by expression outside a spherical region of radius 
Tc- The F core is contained in the region. The boundary condition on the radial components 
ipQ and ipi aX r = couples the two channels. This approach essentially relies on the use of 
spherical coordinates. Here we look for a coordinate independent description of the channel 
mixing that can be efficiently implemented in numerical calculations of the time evolution of 
the system (wave packet propagation). The origin of the kinematic coupling of the channels 
is now sought in symmetries of the system. It is believed that such an approach is rather 
general and could be applied to other many-electron systems where a conventional mean 
field approach does not provide a good approximation. 

The basic idea is that the Hilbert space spanned by isovectors is too large and contains 
states which are physically not acceptable. It is clearly seen already from the fact that the 
radial components of the free ion state @ are not independent. We shall then constrain the 
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state to allow only one bound ion F~ of the form (j^). Any other state with components 
TpQ ~ p(r)Yio(6', yj) and ipi ~ p{r)Yii{6, ip) should be forbidden. This corresponds to making 
the state (H) orthogonal to the vector 



{r\Q) 



( Po{z,p) ^ 



(7) 



Since the state cannot also occur as a virtual (or intermediate) state of the ion in the 
time evolution of the system, we demand that the time dependent wave packet (Q) must be 
orthogonal to the vector \Q): 



{Q\<ij{t)) = {po\Mt))-{Pi\Mt)) = o 



(8) 



for any t > 0, where (po,i|'^o.i(^)) stands for a standard scalar product (written in cylindrical 
coordinates because the functions po,i and ipo^i depend on z and p only). 

From the physical point of view the constraint (§) simply means that there is an unwanted 
scattering mode in our effective one-electron problem. Although the Hamiltonian may allow 
for such a mode, we have given physical reasons to forbid it. The constraint implies 
that the time evolution of the two components of the isovector (^) is no longer independent 
even though the Hamiltonian (|^) does not provide any direct coupling of them. The link 
between the detachment channels with m = and |m| = 1 steams directly from the free ion 
structure. It implicitly assumes that the correlation between the wave functions ipo and ipi 
in the ion perturbed by an interaction with a metal surface is the same as in the free ion. 
Thus, the electronic structure of F~ has been modeled by two one-electron problems linked 
by the constraint (H). 



B. Physical Hamiltonian 

Now we face a problem of incorporating the constraint into quantum dynamics gener- 
ated by the Hamiltonian (^. The difficulty is clear. Suppose an initial state satisfies the 
constraint (||). Applying the evolution operator exp(— itH) to it, we immediately observe 
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that the evolved state fails to satisfy the constraint. The procedure we propose is based on 



the projection operator formalism first introduced for gauge theories p2| , p2| . It has been 
generalized to general constrained systems p3| , |3^ (see also the review [^). The key steps 
are as follows. 

Consider the projection operator 

p = i-|g)(Q|, (9) 

where I is the unit operator, I = for any It is easy to convince oneself that the 
operator is self-adjoint, P"!" = P, and satisfies the characteristic property of a projection 
operator, P^ = P. By construction, the state P |\I') satisfies the constraint @ for any state 
The operator (P) projects any state to the physical subspace defined by the condition 
In particular, P \Q) = 0. To eliminate the state \Q) as a possible intermediate state of 
the system in the time evolution, the Hamiltonian is projected onto the physical subspace 

H ^ PHP = Hphys . (10) 

The physical Hamiltonian (|T0|) is self-adjoint. Hence the time evolution generated by it is 
unitary. The state PHP satisfies the constraint (|^). Clearly, the physical Hamiltonian 
is nonlocal, in general, even if the original Hamiltonian has a standard form of the sum of 
potential and kinetic energies. However, classical limits of H and Hphys are the same. 
The evolution operator has the form 

U(ti,t2) = PTexp (^-i rfrPHP^ P = P exp (-itPHP) P , (11) 

where t = t2 — ti and T exp stands for the time ordered exponential. The second equality 
holds when the Hamiltonian H does not explicitly depend on time. The projection operators 
before and after the exponential in (|1T]) can be omitted if the initial state satisfies the 
constraint (|). 

The physical Hamiltonian provides the sought-for channel mixing. To find terms in the 
new Hamiltonian which give rise to the channel mixing, we compute the action of PHP on 
a generic state A straightforward computation leads to the following result 
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PHP|^) 



Ai 
A2 



^//|^i)-Ai(if + A2) |pi)-A2 \pi) 

(QI^) = (P0|^l) + (Pl|^2), 



(12) 



(QlVsl^^) = {po\Vs\^o) + {pi\Vs\^i) . 



If the state belongs to the physical subspace then, according to (||), Ai = 0. However, 
the action of PHP even on the physical states is not reduced to that of H because the 
amplitude A2 is generally not zero. Thus, after the projection (|10|) the original Hamiltonian 
(I) acquires an additional term 



PHP = (H + W) P 



(13) 



W = -\Q){Q\ Vs. 



(14) 



The operator W provides the channel mixing even if the initial state is in the physical 
subspace. The correlation between the two components of |\1') is a consequence of the 
unbalanced action of the operator Vs on the components of |\E'). It vanishes in the limit of 
the free negative ion. 



C. Interaction potentials 

The potential V^t in the Hamiltonian (P) represents the interaction between the active 
electron and the fluorine neutral core. It is taken as a local model potential which depends 
on the distance r between the electron and the atom center. The potential also includes 
a long range polarization interaction. Its explicit form has been adjusted to reproduce the 
binding energy of the ion F~ as well as the mean radius of the p-orbital. The potential reads 



Kt = -Uo + gr^- ae"^' 



at 



a 

'2^ 



ae 



r < 1 



r > 1 , 



(15) 
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where Uq = 5.64, g = 3.76, a = 1.2558 and the atomic polarizabihty of a fluorine atom 
a = 3.76 All constants are given in the atomic units. 

The potential Vs in the Hamiltonian describes the interaction of the active electron 
with the Cm(111) and Ag{lll) surfaces. It has been proposed by Chulkov et al on the basis 
of their ab initio studies . This local potential depends only on the electron coordinate 
z along the surface normal. Its explicit form can be found in Ref. I^Bf. Qualitatively, 



it is an image charge potential in vacuum which joins smoothly an oscillating potential 
with the period being that of the (111) planes inside the metal bulk. When describing 
an electron motion in the direction perpendicular to the Cm(111) or Ag{lll) surface, this 
potential represents rather well important features of the surface such as the projected band 
gap (between -5.83 and -0.69 eV with respect to vacuum for Cu and between -4.96 and 
-0.66 eV for Ag), the surface state (5.33 eV and 4.625 eV below vacuum for C'u(lll) and 
Ag{lll), respectively), and the image state energy positions (0.82 eV and 0.77 eV below 
vacuum, respectively, for Cu(lll) and Ag (111)). Since the RCT process mainly favors 
transitions around the surface normal, this potential is expected to account for the effect of 
the pecularities of the Cm(111) and Ag{lll) surfaces. 

In order to illustrate the effects of the projected band gap of the C'u(lll) and Ag{lll) 
surfaces, we also present results obtained for a free-electron description of the metal sur- 
face. In this case, the local electron- surface interaction potential corresponds to the y4Z(lll) 
surface and is taken from the work of Jennings et al . 



D. Wave packet propagation 



In the wave packet propagation (WPP) approach |8|,p3|, one studies the time evolution of 



an electron wave function (r |\E'(t)) generated by the system Hamiltonian. Here we consider 
both the static and dynamic problems. In the former, the distance between the ion and the 
metal surface is flxed, while in the latter the ion collides with the surface. In both cases, 
the initial state is the free ion state (^. 
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The time evolution can be regarded as a sequence of infinitesimal time steps generated 
by the evolution operator (111]) 



+ At)) = U(t,t + At) (16) 

Since the initial state |\l/(t = 0)) = 1^1/^-) is orthogonal to the vector \Q), i.e., it is in 
the physical subspace, we can omit the projection operators to the left and right of the 
exponential in (pl^. So in ([TT|) we take 



U(t, t + At)= e-^^PHP _ g-.At{H.,-i?.|Q>(Q|+pvsP) ^ (17) 

where the property Hat \Q) = Ei,\Q) of the vector \Q) has been used; £'a is the eigenvalue 
of i^at corresponding to the eigenf unctions p{r)Yi^m{0,f) with / = 1 and m = 0, ±1. For 
an infinitesimal time step At, the action of the evolution operator (^) can be evaluated by 



means of the split operator approximation p8| , p9 



U(t, t + At) = e-*f PVsPg-^A*(H..-i?.|Q>(g|)g-i^pvsP ^ ^(^^3) _ (^g) 

Making use of the commutation relation of Ha and \ Q) {Q\, this representation can be further 
simplified 

U(t,t + At) = e-*^PVsP ll^g^Aii?. _ _^ j|g-^AtH..g-.AipVsP ^ ^(^^3) (^g) 

The action of the exponential involving Vs is evaluated via a Taylor expansion in which four 
terms are typically kept for the time step At = 0.025 atomic units. 

The action of the kinetic energy operator is computed in the cylindrical coordinates 
which are well suited to the symmetry of the problem 

where Tz contains only the 2;-derivative. The exponential of Hat in ( P^ is then transformed 

as 

g-iAtHat ^ ^-iMH^t J ^ g-i^(T,+yat)g-iAiTpg-i^(T,+yat) J _^ O(At^) . (21) 
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Recall that the operator Hat is diagonal in the isotopic two-dimensional space. So, the 
operator ( [2T| ) acts on both components of |\l'(t)) in the same way. Finally, all the exponentials 
in (|2T| ) are approximated by means of the Cayley representation pO[ 



-lAtA 



1 



■At 



A 



(22) 



In order to accurately reproduce the wave packet variation close to the atom center, we use 
a mapping procedure |4l|,^ defined by 

0.95^3 



z = f{0 =0.05^ + 

p = fin) = 0.057] + 



400 + ^2 
0.95r/3 
400 + 7] 



(23) 



2 ' 



The wave packet ^{t,z,p) is evaluated on a 2D mesh of points {C,k,Vj) of the size 
1200x800 with the step size A = 0.2 atomic units for both the coordinates. At the grid 
boundary, an absorbing potential is introduced in order to avoid the wave packet 

reflection. 

The kinetic energy operator ( pO|) has to be written in the auxiliary variables ^ and t] and 



then discretized. After the change of variables (|23D the operator Tp assumes the form 

1 1 d f d 1 



1 



(24) 



2J,/fdvJdvVf 2/2 
where J{ri) = f'{rj) is the Jacobian. The grid in the //-coordinate is set as rjj = A/2+A(j — 1). 
For every value ^ we have = \l/(^, r/j), and the action of ( p4D is defined by the following 
midpoint procedure 



2A2 j.^/T. 



fj+1/2 



j+1 



X 



i+i/2 



J. 



j-l/2 



+ 



1 

2lf 



(25) 



Here the subscript j ± 1/2 means that the corresponding function is taken at the midpoint 
?7j ± A/2. A similar expression can be obtained for the action of on the grid C,k = 
eo + A(A;-l). 
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In the first series of calculations, we study the static problem when the ion F~ is at a fixed 
distance Z from the metal surface. The survival amplitude of the ion (the auto-correlation 
of the wave function) is given by 



From the Laplace transform of the function A{t), one can obtain the density of states (DOS) 
projected on the free ion wave function. The structure of the DOS yields the energy level and 
its width for the negative ion state interaction with the surface. In what follows, this width 
is referred to as the "static width" to emphasize that it is extracted from static calculations. 
It gives the electron transfer rate between the negative ion and the metal surface in the 
static problem. 

In the second series of calculations, we study the evolution of the electron wave packet 
when a negative fiuorine ion collides with the surface. The ion is assumed to approach the 
surface along a straight line perpendicular to the surface at a constant velocity v. Only 
the incoming part of the collision is studied. The time dependence of the wave function 
is obtained in the projectile reference frame, i.e., the time dependence of the Hamiltonian 
occurs through the potential Vs. For each collision velocity, the ion survival probability, 
P{t,v) = |74(t, v)p, is computed. To analyze the dynamics of the charge transfer, we define 
an effective width of the negative ion state by 



where Z — Zq — vt with Zq being an initial distance of the ion from the metal surface. 
It corresponds to an effective decay rate of the ion when it approaches the surface with a 
velocity v. Comparing G{Z,v) to the level width obtained in the static calculations allows 
us to see to what extent the dynamical evolution can be described by the static width of 
the ion level. 




(26) 





(27) 



dt 
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III. RESULTS AND DISCUSSION 



A. F ions interacting with a surface Al(lll) 

The interaction of an F" ion with an ^4/(111) surface, where the latter is regarded as a 
free-electron metal surface, has already been studied by the CAM method associated with 
the effective range treatment of the negative ion It lead to a successful description 



of the negative ion formation in a grazing angle scattering |]I2|. Similarly, for large angle 
scattering from Ag{llQ), and polycrystalline Ag and Al surfaces, a quantitative agreement 
with experiment results [p!l| , ^ has been obtained [Q. 



In Figures 1 and 2 we compare the results obtained in the static case by two different 
methods: The CAM method with effective range treatment of the negative ion and the 
present WPP results obtained with the projection formalism. In both cases, the AZ(lll) 
surface is described as a free-electron metal surface using the potential proposed in . 



Figures 1 and 2 present, respectively, the energy position and the width of the F~ ion level 
interacting with the surface as a function of the ion-surface distance Z measured from the 
image plane. The characteristics of the ion level as a function of Z display the behavior 
common for atomic species in front of a free-electron surface: The energy of the negative ion 
state decreases as the ion is placed closer to the surface, which can be anticipated because of 
the image charge attraction; The level width increases roughly exponentially as Z decreases. 
The results obtained by two different methods are extremely close to each other. This gives 
confidence in the equivalence of the two descriptions of the F~ ion. The results for a free- 
electron y4/(lll) surface are used below as a "free electron" reference to which compare the 
C'u(lll) and Ag{lll) results. It appears that the free electron results are almost identical 
for the three metals, except at very small distances from the surface. 
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B. F ions interacting with a surface Cu(lll). A static case 

Figures 3 and 4 present the F~ ion level characteristics (energy (Fig.3) and width (Fig. 4)) 
as a function of the ion-CM(lll) surface distance. The negative ion level energy exhibits 
an avoided crossing around 4ao from the surface, which is quite different from the smooth 
behaviour seen in Figures 1 and 2 for the free-electron metal. This is a direct consequence 
of the peculiarities of the Cm(111) surface and a similar situation has already been observed 
in the case of H- interacting with the same surface ||23|| . 

A schematic picture of the electronic structure of the model Cm(111) surface is shown in 
Fig. 5. The energy of electronic levels is plotted as a function of the electron momentum, fcy, 
parallel to the surface. For k\\ equal to zero, the projected band gap lies within the energy 
range from -5.83 to -0.69 eV (with respect to vacuum). Inside the gap, there is a surface 
state at -5.33 eV. In the present model of a Cu surface, the dispersion curves of all the metal 
electronic states as functions of k\\ are parabolic with a free-electron mass. 

The resonant charge transfer process corresponds to transitions between the ion level 
and metal states of the same energy. At large distances, the F~ ion level is degenerate 
with the band gap. Therefore it can only decay to metal states with a finite k\\ i.e. either 
to the 2D surface state band, or into 3D propagating states. As Z decreases, the energy 
of the negative ion state decreases and it comes close to the bottom of the 2D surface 
state continuum. The ion can decay by ejecting an electron with the angular momentum 
m = 0, ±1 (the quantization axis is normal to the surface). As explained in Ref. [^, a 



resonance cannot cross the bottom of a 2D continuum in the symmetrical case m = and 
there always exists a bound state below the bottom of the continuum. This state has an 
avoided crossing with the state which becomes the free ion state as Z tends to infinity. The 
F~ ion character is then found to be associated with two different states depending on the 
Z range. It is transferred from the upper to the lower state when going through the crossing 
region (decreasing Z). Far from the avoided crossing region, the ion energy level is rather 
close to that found in the free electron case. 
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As for the width, at large Z, its absolute value for F~ in front of the Cu{lll) surface is 
larger than that in the case of a free-electron surface. This result might appear surprising 
since the projected band gap prohibits the electron transfer from the projectile to the metal 
along the surface normal (fcy = 0) and blocks the RCT into the 3D bulk continuum. Indeed, 
as can be seen in Fig. 5, there are no electronic states of the metal with small k\\ which 
are in resonance with the negative ion state. The potential barrier separating the ion and 
the surface attains its least value in the direction normal to the surface. Therefore, the 
surface normal is the preferred direction of the resonant electron transfer. One would then 
expect that the effect of the projected band gap would be to stabilize the negative ion level 
as compared with the case of a free electron metal. This has indeed been found for H~ 
interacting with C'u(lll) where the width of the H~ state was much reduced as compared 
to the /Al{lll) - case [|3 . 



In contrast, we have observed an increase of the fluorine negative ion decay rate as 
compared to the free-electron metal case. The reason is twofold. First, the 2D surface state 
continuum contributes to the decay of F~ . Second, a fluorine has a much larger electron 
affinity than a hydrogen. Thanks to a better overlap of the wave functions, the 2D surface 
state continuum, when energetically allowed, is a dominating decay channel for an ion state 
lying within the band gap p3ip8| , ^ . Moreover, when the binding energy of a negative ion 
is close to that of the surface state, a coupling of the ion level with the 2D surface state 
continuum is more efficient than its coupling with the 3D continuum of the free-electron 
metal states (see a discussion in [^]). The efficiency of the 2D surface state continuum as a 
decay channel can also be deduced from the sharp decrease of the level width when passing 
through the crossing region as Z decreases, i.e., when this decay channel becomes closed. 

For small values of Z: (i) the 2D surface state continuum does not contribute to the 
decay of the negative ion, and (ii) the energy of the negative ion state is very close to the 
bottom of the projected band gap for small k\\ so that the band structure effect for the decay 
into the 3D bulk continuum vanishes. Therefore, we find the width of the level with ionic 
character very close to the free-electron results. 



Finally, we can stress that the procedure of extracting the resonance characteristics 



employed here is based on the autocorrelation function p6| ) using the free F~ ion wave 
function as the initial state. It converges well for the states of an ionic type. Convergence 
of the resonance characteristics for other states is difficult to achieve. This is the reason for 
showing only one state far from the crossing region (small or large Z) . In the crossing region, 
the ionic character is shared between the lower and upper states so that the characteristics 
(energy and width) of both of them can be extracted. Since the convergence is easier to 
achieve for the energy of the state, the interval of distances Z where both states are presented 
is larger in Fig. 3. 



C. F ions interacting with an Ag(lll) surface. Static and dynamic studies 

The electronic structures of Ag{lll) and Cu(lll) look rather similar (cf. Figures 5 and 
6). However, characteristic features of the electronic structures occur at different energies. 
The surface state in Ag{lll) is located higher in energy than in Cm(111). For this reason 
the avoided crossing appears at a larger Z where the ion-surface charge transfer interaction 
is smaller. As a consequence the avoided crossing could not be resolved in the energy 
dependence because it is localised in a too small range of Z. Therefore, we have chosen to 
represent the results for the energy and the width by a single continuous line (Figures 7 
and 8, respectively). In fact, the energy of the negative ion state is almost the same as for 
the free-electron metal surface. The characteristic change of the level width when passing 
the crossing region (decreasing Z) is however still perfectly visible. It fully confirms the 
dominance of the 2D surface state channel in the F~ ion decay at large Z. 

At small Z, the F~ ion level is embedded in the 3D propagating states of Ag and its 
characteristics are very close to those of the free-electron case. Very close to the surface, 
the decrease of the level width as compared to ^4/(111) free-electron results is caused by the 
closeness of the F~ level to the bottom of the Ag{lll) valence band. As can be seen in Fig. 
3b the bottom of the Ag{lll) valence band is located at —9.7 eV in our model description 
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of the Ag{lll) surface. For the model free-electron case it is located at —15.9 eV. 

By studying the time dependent problem, one can find out whether the peculiarities of 
Ag{lll) observed in the static case can still be visible in a collision of F~ with the Ag{lll) 
surface. We have computed the effective level width G, defined by (0), as a function of 
the distance for an ion F~ approaching the surface at different velocities. The results are 
displayed in Fig. 9 together with the results of the static case for Ag{lll) and the free- 
electron metal surface AI{111). As the collision velocity is increased, the effective width 
becomes closer to the free-electron result. This feature is very similar to what we have 



found for ions H~ interacting with a Cm(111) surface ||2^. The system needs a finite time 
to react on the presence of the projected band gap. If the collision is too fast, the electron 
wave packet does not have enough time to "explore" the band structure of the metal, and 
the ion decay remains identical to that on a free-electron metal surface. 

In contrast, as the collision velocity is decreased, the effective width comes nearer to 
the static Ag{lll) width. For the smallest velocity used here, 0.0058 atomic units which 
correspond to a collision energy about 16 eV, the effective width is very close to the static one 
at large Z. When Z is decreased, the effective width fails to perfectly reproduce the change 
of the behavior associated with the crossing of the bottom of the surface state continuum 
for all collision velocities considered here. We can nevertheless see that this variation is 
better reproduced as the collision velocity is decreased. In fact, the dynamical broadening 



introduced by the change of the negative ion state energy with time makes it impossible 



to reproduce a sharp variation of the static width as the ion approaches the surface. This 
leads to rounded delayed structures displayed in Fig. 9. Moreover, the oscillations of the 
effective width at small Z can tentatively be attributed to the population transfer between 
the two adiabatic states at the crossing region. From these results one can conclude that in 
the studied collision energy range, 16 eV - 5 keV, the charge transfer rate is intermediate 
between the free-electron case and the static Ag{lll) case. 

The formation of F~ ions and ions by collision on a Ag (111) surface has been studied 
experimentally by Guillemot and Esaulov . When comparing with results obtained with 
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free-electron like surfaces they found that the H~ data presents a strong band gap effect. 
In particular, the survival probability of the negative ions leaving the surface was much 
larger for Ag{lll). This was attributed to the blocking of the resonant charge transfer in 
H~ /Ag{lll) system, in line with theoretical results obtained for H~ /Cu{lll) [^] where 
the RCT rates are reduced by orders of magnitude compared to the free-electron case. 
At the same time, results obtained with F~ ions were approximately consistent with a 
description based on the charge transfer rate obtained in the framework of a free-electron 
description of the metal surface. In the energy range studied in their work, as shown above, 
the dynamical behaviour of the charge transfer is intermediate between the free-electron and 
static Ag{lll) limits. This prohibits the use of the simple classical treatment of the parallel 
velocity effect [p8| , ^ which was shown to be important for a formation of F~ , even at rather 
low collision energy Consequently, we cannot quantitatively compare our results with 
theirs. Nevertheless, we can notice that in our study the band gap effect is reversed and much 
smaller in the present system than in the H~ /Cu{lll) system; it is even partly suppressed 
by the ion motion. These findings qualitatively agree with the experimental results . 



IV. CONCLUSIONS 

We have reported on a study of the electron transfer in the ion-metal systems 
/Cu{lll) and F^ /Ag{lll). We have developed a new method to describe the effect 
of six quasi-equivalent outer-shell electrons of F~ on the resonant charge transfer process. 
The original six-electron problem has been transformed into two one-electron problems in 
which the dynamics are not independent but rather are linked via a constraint. The projec- 
tion formalism for quantum systems with constraints has been used to obtain the quantum 
mechanical propagator for such a system. This modeling of the ion F~ is simple, efficient 
and can easily be implemented in the wave packet propagation approach. 

Both the Cm(111) and Ag{lll) surfaces exhibit an electronic structure with a projected 
band gap in which the ion energy level lies at large ion-surface distances. This peculiarity 



21 



of the electronic structure influences the charge transfer interaction with the ion, leading to 
a few remarkable features: 

• The ion level presents an avoided crossing with the bottom of the surface state con- 
tinuum as predicted for a 2D continuum with m = 0. The avoided crossing is very clearly 
marked for the system /Cu{lll). 

• Because of the correlation between the six electrons of the F~ ion the avoided crossing, 
which is a characteristic feature of the symmetric case m = 0, appears in the present system 
where electrons in both states with m = and |m| = 1 contribute to the charge transfer. 

• When the negative ion level is low in the projected band gap, the band gap does not 



cause a drastic drop of the charge transfer rate as observed in other systems p3|,^. This 
feature of the charge transfer has been attributed to the following effects: (i) The band gap 
effect is expected to decrease as the projectile level is lower in the band gap; (ii) The decay 
to the surface state is favored over the decay to 3D bulk states because of a greater spatial 
overlap of the electron wave function with the surface state; and {Hi) The polarization of 
negative ions does not enhance the band gap stabilization effect as it does for neutral atoms 
(see a discussion in Ref. P^). 

• At small ion-surface distances, when the negative ion state is not in the band gap or 
inside the gap but close to its bottom, the energy and the width of the negative ion level 
are practically identical to those found in the free-electron surface case. 

• Studies of the corresponding dynamical systems, when the ion F~ approaches the 
surface, have shown that the above features survive, although partially, over a large collision 
energy range. 
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FIG. 1. Energy position of the F~ ion level in front of the AZ(lll) free-electron-like surface, 
as functions of the ion-surface distance, measured from the image plane (atomic units). The solid 
line represents the results obtained with the present WPP approach. The black dots indicate the 
results obtained with the CAM method. 
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FIG. 2. Energy width for the same system as in Fig.l 
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FIG. 3. Energy position of the F~ ion level in front of the model Cu{lll) surface, as functions 
of the ion-surface distance, measured from the image plane, (atomic units) The energy reference 
is the vacuum level. Black dots: results for the free-electron Al(lll) surface. The horizontal 
dashed-dotted line indicates the energy position of the bottom of the surface state continuum. Solid 
line: results for the highest lying resonance. Dashed line: results for the lowest lying resonance. 
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FIG. 4. Energy width for the same system as in Fig.3. 



29 




FIG. 5. A schematic picture of the electronic structure of Cu(lll) (work function 4.9 eV) as 
a function of the electron momentum parallel to the surface (atomic units). The energy reference 
is the vacuum level. The shaded area represents the 3D valence band continuum. The dashed 
line represents the 2D surface state continuum. The energy of F~ level at some distance from the 
surface is displayed as the horizontal solid line. 
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FIG. 6. Same as Fig. 5 for the ^^(111) surface (work function 4.56 eV). 
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FIG. 7. Energy position of the F~ ion level in front of the model ^^'(111) surface (solid lines), 
as functions of the ion-surface distance, measured from the image plane (atomic units) . Black dots 
represent the results obtained for the free-electron Al(lll) surface. The horizontal dashed-dotted 
line indicates the energy position of the bottom of the surface state continuum. 
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FIG. 8. Energy width for the same system as in Fig. 7. 
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FIG. 9. The effective level width G versus the ion-surface distance Z for various collision 
velocities. The solid dots and solid squares stand for, respectively, the free electron and ^^(111) 

static widths. Continuous lines represent the results obtained by the constrained wave packet 
propagation method for the model surface ^^'(111) for various collision velocities (see insert). 
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